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Phosphatidylcholine  Monolayer  Formation  at  a  Liquid:  Liquid  Interface 
as  Monitored  by  the  Dynamic  Surface  Tension 

Robert  A.  Walker  and  Geraldine  L.  Richmond* 

University  of  Oregon,  Eugene,  OR  97403 

Abstract:  Dynamic  surface  tension  experiments  have  monitored  the  rate  of  phospholipid 
monolayer  formation  at  a  planar  aqueous:CCl4  interface  from  a  solution  of 
phosphatidylcholine  (PC)  vesicles.  The  rates  at  which  monolayers  form  show  that 
monolayer  formation  is  a  barrier  controlled  rather  than  a  diffusion  controlled  process.  At 
sufficiently  low  bulk  PC  concentrations,  monolayer  formation  kinetics  are  first  order  in 
aqueous  PC  concentration.  The  kinetics  can  be  described  by  a  rupture  mechanism  which 
postulates  that  vesicles  disintegrate  at  the  aqueous:CCl4  boundary  allowing  all  of  the 
monomers  to  adsorb  to  the  interface.  When  coupled  to  an  appropriate  two  dimensional 
equation  of  state  this  mechanism  quantitatively  models  the  observed  dynamic  surface 
tension  data.  Solutions  of  vesicles  in  their  liquid  crystalline  phase  form  tightly  packed 

monolayers  at  concentrations  above  ~2  pMolar,  while  solutions  of  gel  phase  vesicles  form 

expanded  monolayers  regardless  of  bulk  PC  concentration.  From  the  Arrhenius  behavior 
of  the  dynamic  surface  tension  data,  we  find  that  the  barrier  to  monolayer  formation  is 
larger  for  solutions  of  liquid  crystalline  vesicles  than  for  solutions  of  gel  state  vesicles. 
These  phenomena  -  the  extent  of  monolayer  formation  and  Ea  for  monolayer  formation  - 
support  a  model  which  describes  monolayer  formation  as  a  thermodynamically  driven 
transformation  from  monomers  within  a  bilayer  to  monomers  adsorbed  at  the  aqueous:CCl4 


interface. 


I.  Introduction 


A  distinguishing  characteristic  of  long  chain  phospholipids  is  their  ability  to  form 
vesicles  in  aqueous  solutions.  These  uni  or  multilamellar  bilayer  structures  serve  as  the 
simplest  biomimetic  models  of  biological  cell  systems  and  considerable  effort  has  gone  into 

studying  how  bilayer  composition  affects  physical  properties  such  as  permeability, *»2 
structural  disorder,3'6  and  chemical  stability.7  Phospholipid  vesicles  have  also  found 
application  as  vehicles  in  drug  delivery  systems.8-9  Given  the  exquisite  control  researchers 

can  exercise  over  vesicle  size  characteristics,10-12  phospholipid  vesicles  appear  to  provide  a 
very  attractive  means  of  transporting  pharmaceuticals  in  vitro  to  specific  areas  within  the 
body.  Essential  to  their  role  as  biological  membrane  models  and  microscopic  transport 
systems  is  the  issue  of  vesicle  structural  integrity.  For  almost  two  decades  researchers 
have  known  that  solutions  of  phospholipid  vesicles  form  monolayers  of  phospholipid 

monomers  at  ainaqueous  interfaces.13-14  Efforts  to  structurally  modify  vesicles  by 
incorporating  membrane  proteins  and/or  cholesterol  have  led  to  formation  of  binary 

monolayer  systems  at  a  variety  of  interfaces.15'17  Related  work  has  examined  intervesicle 

exchange  of  phospholipid  monomers.18'24  These  phenomena  -  monolayer  formation  and 
monomer  exchange  -  necessarily  require  that  vesicles  themselves  be  dynamic  structures 

rather  than  rigid,  static  bodies.25  The  goal  of  this  work  lies  in  examining  phospholipid 
monolayer  formation  at  a  planar  liquiddiquid  interface  where  one  of  the  fluid  phases  is  an 
aqueous  solution  of  phosphatidylcholine  (PC)  vesicles.  With  the  Wilhelmy  plate  method 
we  monitor  the  rate  of  monolayer  formation  by  measuring  the  dynamic  surface  tension  of 
an  aqueous:carbon  tetrachloride  interface.  We  show  that  monolayer  formation  kinetics 
depend  not  only  on  aqueous  PC  concentration  but  also  on  vesicle  bilayer  phase.  We  also 
calculate  activation  energies  for  the  monolayer  formation  process  and  find  that  vesicles  in 
the  liquid  crystalline  phase  must  surmount  a  larger  barrier  to  form  monolayers  than  gel 
phase  vesicles. 
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In  1978  Verger  and  coworkers  demonstrated  that  monolayers  of  phospholipid 

monomers  will  form  from  aqueous  solutions  of  phospholipid  vesicles.13  Shortly  thereafter 
Schindler  examined  monolayers  formed  in  this  fashion  at  the  aqueous:air  interface  and 

found  that  monomers  exchange  freely  with  an  underlying  layer  of  phospholipid  vesicles.14 
These  studies  concluded  that  the  surface  layer  possesses  an  appreciably  higher 
concentration  of  vesicles  than  the  bulk  solution  and  that  in  the  limit  of  an  infinitely  dilute 

monolayer,  vesicles  probably  disintegrate  completely  on  contact  with  the  interface.14 
MacDonald  and  Simon  examined  the  relationship  between  lipid  bilayer  phase  and  the  rate  of 

monolayer  formation.26  Working  with  solutions  of  DMPC  vesicles  above  the  gel-liquid 
crystalline  transition  temperature  (Tc  =  23°  C),  the  authors  found  that  aqueous  solutions  of 
DMPC  vesicles  rapidly  equilibrated  with  the  airiaqueous  interface,  forming  monolayers 
with  a  surface  pressure  of  49  mN/m.  In  contrast  aqueous  solutions  of  DMPC  vesicles 
below  Tc  required  days  to  equilibrate  and  even  then  the  monolayer  formed  was 
considerably  more  expanded  than  that  formed  from  vesicles  above  Tc. 

Similar  effects  have  been  observed  in  rates  of  intervesicle  monomer  exchange. 
McLean  and  Phillips  noted  that  the  rate  of  DMPC  monomer  exchange  between  small, 
unilamellar  DMPC  vesicles  exhibited  an  abrupt,  discontinuous  rise  in  the  vicinity  of  the 

lipid  bilayer  gel-liquid  crystalline  transition  temperature.22  Additional  studies  reported  that 
monomer  exchange  between  vesicles  constituted  a  first  order  process20  and  the  rate  of 

transfer  decreased  exponentially  with  increasing  acyl  chain  length  27  Two  different 
mechanisms  arose  to  explain  these  results:  a)  transfer  of  monomers  through  solution  and  b) 
monomer  transfer  via  vesicle  collision.  Yang  and  Huestis  demonstrated  how  each 
mechanism  provided  a  consistent  explanation  of  observed  transfer  phenomena  under 

different  conditions.18  Through-solution  exchange  adequately  described  data  for  more 
hydrophilic  lipids  while  the  collision  mechanism  worked  well  for  more  hydrophobic 
species.  Potentially,  both  mechanisms  could  play  important  roles  during  formation  of  a 
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monolayer  from  a  solution  of  PC  vesicles.  If  monolayer  formation  proceeded  primarily  by 
means  of  through-solution  transport  of  monomers,  then  we  might  expect  monolayers  to 
form  quite  slowly  as  monomers  adsorbed  to  the  interface  one  by  one.  However,  a 
mechanism  which  required  a  collision  between  vesicles  and  the  interface  would  favor  the 

vesicle  disintegration  picture  inferred  in  earlier  studies14  and  allow  monolayers  to  form 
more  rapidly. 

The  experiments  described  in  this  paper  address  the  way  in  which  aqueous 
solutions  of  phosphatidylcholine  vesicles  form  monolayers  of  monomers  adsorbed  to  a 
planar  aqueousxarbon  tetrachloride  interface  (Figure  1).  Specifically,  we  Find  that  the  rate 
of  monolayer  formation  depends  on  bulk  PC  concentration  as  well  as  temperature. 

Dynamic  surface  tension  experiments  show  several  different  regimes  of  monolayer 

formation  kinetics.  At  bulk  PC  concentrations  below  ~2|iMolar  and  at  temperatures  above 

Tc ,  the  surface  tension  exhibits  single  exponential  decay  behavior,  slowly  approaching  an 
asymptotic  limit  characteristic  of  the  bulk  PC  concentration  and  the  system’s  temperature. 

In  the  high  concentration  limit  above  Tc,  the  surface  tension  diminishes  rapidly  in  a  fashion 
not  readily  characterized  by  a  simple  functional  form.  When  the  experimental  temperature 
lies  below  a  system’s  Tc,  monolayer  formation  happens  very  slowly,  regardless  of  bulk  PC 
concentration.  We  show  how  a  simple  kinetic  model  coupled  to  an  appropriate  two 
dimensional  equation  of  state  explains  all  of  the  observed  phenomena.  The  mechanism 
requires  that  PC  monolayers  form  from  complete  rupture  of  vesicles  at  the  interface.  A 
mechanism  which  postulates  that  monomers  adsorb  individually  does  not  accurately 
describe  the  observed  kinetics.  From  studying  the  temperature  dependence  of  the 
monolayer  formation  kinetics,  we  calculate  activation  energies  for  the  monolayer  formation 
process  and  find  fundamentally  different  behavior  depending  on  whether  aqueous  phase 
vesicles  are  in  their  gel  or  liquid  crystalline  phase.  Our  results  share  similarities  with  data 
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from  intervesicle  monomer  exchange  studies,  and  we  can  relate  our  findings  to  physical 
properties  of  the  vesicle  bilayer  structure. 


II.  Experimental 
A.  Sample  Preparation 

The  phosphatidylcholines  used  in  this  work  belong  to  a  family  of  saturated, 
symmetric,  diacyl  phosphatidylcholines  (Figure  2a)  with  chains  ranging  in  length  from  C]2 
(dilauroyl-sn-phosphatidylcholine,  DLPC)  to  C18  (distearoyl-sn-phosphatidylcholine,  .. 
DSPC)  in  increments  of  two  methylene  units  (C14  =  dimyristoyl-sn-phosphatidylcholine, 
DMPC;  C16  =  dipalmitoyl-sn-phosphatidylcholine,  DPPC).  In  aqueous  solution  these 
phosphatidylcholines  spontaneously  aggregate  to  form  bilayer  sheets.  Upon  sonication, 
these  sheets  reorganize  to  form  uni  or  multilamellar  bilayer  vesicles.12  The  thermodynamic 
phase  of  these  bilayer  vesicles  depends  on  the  temperature  of  the  solution:  at  low 
temperatures  the  bilayer  interiors  are  frozen  in  a  gel  phase  (Pp)  while  at  higher 
temperatures,  the  chains  melt  and  assume  a  more  fluid,  disordered  liquid  crystalline  phase 
(LJ.  (Figure  2b)  Phase  transition  temperatures  (T^  vary  with  PC  chain  length  (Table  1). 

The  solubility  of  PC  monomers  in  solution  remains  exceedingly  small  during  all  of 
the  experiments  described  below.  Using  fluorescent  tagged  PC  monomers,  Nichols 
measured  the  bilayenfree  monomer  partition  coefficients  for  analogs  of  DMPC,  DPPC  and 
DSPC  and  found  partitioning  into  the  bilayer  to  increase  sharply  with  each  additional  pair 
of  methylene  units  (~1  x  107  for  DMPC  analog,  ~1  x  108  for  DPPC  analog,  ^  109  for 
DSPC  analog). 21  For  our  experiments  aqueous  PC  concentration  lies  between  ~0.5  and 
10  jiMolar,  meaning  that  even  at  the  highest  aqueous  concentrations,  the  number  of  free 
DLPC  monomers  in  solution  represents  less  than  1%  of  the  number  necessary  to  form  a 
tightly  packed  monolayer  (extrapolating  the  partition  coefficient  for.DLPC  to  be  1  x  106 
bilayer/free  monomer). 
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The  four  different  phosphatidylcholines  were  purchased  as  lypholized  powders 
from  Avanti  Polar  Lipids  (>99.9%  pure)  and  used  as  received.  The  aqueous  phase 
consisted  of  H20  (Bamsted  Nanopure,  18  M£2)  that  had  been  washed  with  CC14  to  remove 
any  small  organic  contaminants.  All  aqueous  solutions  were  buffered  with  a  sodium 
phosphate  solution  to  a  pH  of  7.0  (10  mMolar  in  phosphate  ion).  The  purity  of  carbon 
tetrachloride  (Aldrich,  reagent  grade  or  Baker,  reagent  grade)  was  confirmed  by  FTIR  and 
'H  NMR.  Interfacial  tension  measurements  of  the  neat  watercarbon  tetrachloride  interface 
produced  results  in  agreement  with  literature  values, 28  indicating  an  absence  of  any 
detectable  surface  active  contaminants  in  the  two  solvents. 

Phosphatidylcholine  stock  solutions  were  prepared  by  suspending  a  measured 
amount  of  phosphatidylcholine  in  5  ml  of  buffered  H20  and  then  sonicating  the  solution 
above  the  appropriate  phase  transition  temperature.  This  procedure  created  a  solution  of 
vesicles  ranging  in  diameter  from  60  -150  nm  (sample  dependent)  as  determined  from 
dynamic  light  scattering  measurements.  Typical  concentrations  of  these  stock  solutions 
ranged  from  0.5  mMolar  to  2.0  mMolar.  Phosphatidylcholine  sample  solutions  used  in  the 
experiments  were  prepared  by  adding  a  small  volume  of  the  stock  solution  to  a  volumetric 
flask  containing  25  mL  of  buffered  H20.  Concentrations  of  the  solutions  used  in 
experiments  varied  between  0.05  (jMolar  and  50  pMolar.  Flasks  used  for  sample 
solutions  were  rinsed  repeatedly  with  PC  solutions  of  equivalent  concentration  in  order  to 
eliminate  any  possible  effect  of  PC  adsorption  to  flask  walls. 

Solutions  prepared  from  freshly  sonicated  phosphatidylcholine  stock  solutions 
showed  higher  surface  activity  than  those  made  from  stock  solutions  which  had  aged 
several  hours  in  agreement  with  the  findings  of  Qiu  and  MacDonald.2^  After  creating  a 
phosphatidylcholine  stock  solution,  we  allowed  a  minimum  of  10-12  hours  to  elapse  before 
preparing  sample  solutions.  Further  waiting  led  to  no  additional  change  in  PC  surface 
activity.  Stock  solutions  were  replaced  every  5  days. 
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B.  Interfacial  Tension  Measurements 

Interfacial  tension  measurements  at  the  watencarbon  tetrachloride  interface  were 
carried  out  with  a  Wilhelmy  plate  microbalance  assembly  (KSV  Instruments)  and  a  6  cm 
diameter  glass  crystallizing  dish.  An  experiment  consisted  of  first  depositing  a  20  ml 
aliquot  of  CC14  into  the  dish  followed  by  a  gentle  25  ml  addition  of  a  PC  solution.  With  the 
interface  established,  a  platinum  plate  was  lowered  so  that  its  top  edge  lay  submerged 
beneath  the  airraqueous  surface  and  its  bottom  edge  extended  ~1  mm  below  the  nominal 
aqueous:CCl4  interface  (Fig  3a).  This  experimental  approach  had  the  disadvantage  of  .. 
missing  the  initial  moments  (-45-60  sec)  of  monolayer  formation.  However,  monolayers 
typically  required  many  hours  to  equilibrate  and  for  sufficiently  low  PC  concentrations,  the 
interfacial  tension  at  the  start  of  data  collection  differed  from  the  tension  of  the  neat 
aqueous:CCl4  interface  by  less  than  1  mN/m.  This  technique  of  premixing  the  PC 
solutions  reduced  the  effects  of  convection  and  anomalous  concentration  gradients  in  the 
aqueous  phase.  We  also  conducted  a  series  of  experiments  in  which  we  first  established  a 
neat  interface  and  then  made  a  small  volume  addition  of  the  high  concentration  PC  stock 
solution,  raising  the  aqueous  concentration  to  the  desired  value.  30-32  This  method  allowed 
us  to  ascertain  the  tension  of  the  neat  interface  thus  ensuring  interfacial  purity.  However, 
this  analytical  convenience  came  at  the  expense  of  a  uniformly  mixed  aqueous  phase 
following  the  small  volume  stock  solution  addition.  In  the  low  concentration  limit  data 
from  both  techniques  agreed  quite  well.  At  higher  bulk  PC  concentrations  the  small 
volume  addition  approach  led  to  an  initial  steep  drop  in  tension  followed  by  a  more  gradual 
approach  to  an  asymptotic  limit. 

Most  experiments  were  carried  out  under  ambient  conditions  (21-23°  C). 

Thermally  controlled  experiments  designed  to  investigate  the  effect  of  temperature  on 
monolayer  formation  kinetics  involved  attaching  a  water  jacket  around  the  crystallizing 
dish.  A  circulating  bath  allowed  for  temperature  control  to  ±  1°  C.  Although  CC14  boils  at 
-70°  C,  the  aqueous  layer  effectively  reduced  the  boiling  point  to  -50°  C  imposing  an 
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upper  limit  to  the  accessible  temperature  range.  Temperature  dependent  dynamic  surface 
tension  measurements  were  carried  out  in  a  window  between  5°  C  and  45°  C. 

As  the  PC  monolayer  formed  at  the  aqueous:CCl4  interface,  the  surface  tension  (y) 
decreased,  slowly  approaching  an  equilibrium  value.  Terminal  y(t)  readings  were  recorded 
when  the  rate  of  change  of  y(t)  fell  to  less  than  2%  of  y0  -  y(t)  per  hour.  (y0  represents  for 
the  tension  of  the  neat  aqueous:CCl4  interface,  44.5  mN/m.28)  Values  of  equilibrium 
surface  tension  were  converted  to  surface  pressures  (ITC  =  y0  -  ye)  and  plotted  versus  bulk 
PC  concentration  in  order  to  map  out  adsorption  isotherms  and  measure  the  extent  of  .. 
monolayer  formation.  Representative  isotherms  for  DLPC  and  DSPC  at  room  temperature 
appear  in  Figure  3b.  Terminal  surface  pressures  of  DLPC  monolayers  adsorbed  to  the 
water:carbon  tetrachloride  interface  agree  reasonably  well  with  reported  results  of 
phosphatidylcholine  monolayers  adsorbed  to  alternative  aqueous:organic  interfaces.33'36 
Analyzing  the  isotherm  data  by  means  of  the  Gibbs  equation  shows  that  aqueous  solutions 
of  DLPC  vesicles  at  room  temperature  form  tightly  packed  monolayers  with  surface 
concentrations  of  1.8  ±  0.3  x  1014  molecules/cm2  (55  A2/molecule).30-31>37  Furthermore, 
this  terminal  surface  coverage  occurs  at  bulk  DLPC  concentrations  of  >  2jxMolar.  In 
contrast  solutions  of  DSPC  vesicles  under  ambient  conditions  form  expanded  monolayers 
regardless  of  bulk  concentration.  Differences  between  the  two  systems  can  be  traced  to  the 
lipid  bilayer  phase  dependent  behavior  of  monolayer  formation.  (Vide  infra.) 

III.  Dynamic  Surface  Tension  Modeling 
A.  Barrier  versus  Diffusion  Controlled  Kinetics 

Solutions  of  phosphatidylcholine  vesicles  form  monolayers  of  phosphatidylcholine 
monomers  at  interfaces.  13-15,26  As  our  primary  objective  we  seek  to  understand 
monolayer  formation  at  the  planar  boundary  formed  between  an  aqueous  solution  of  PC 
vesicles  and  an  underlying  organic  phase  of  carbon  tetrachloride.  The  issue  of  surfactant 
adsorption  kinetics  has  received  considerable  attention  for  many  years.38'45  In  particular  a 
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number  of  studies  have  specifically  examined  adsorption  kinetics  at  liquidiliquid 
interfaces.46'50  Adsorption  can  be  described  as  having  two  limits:  diffusion  controlled 
adsorption  and  barrier  controlled  adsorption.  Both  situations  define  a  subsurface  layer  as 
an  intermediate  region  between  the  surface  and  the  bulk  solution.45’47  In  the  case  of 
diffusion  controlled  kinetics,  adsorption  from  the  subsurface  to  the  surface  occurs  quickly 
and  diffusion  of  surfactants  from  the  bulk  to  the  subsurface  layer  represents  the  limiting 
step  in  monolayer  formation.  The  rate  of  monolayer  formation  obeys  traditional  laws  for 
mass  transport  and  the  subsurface  concentration  remains  depleted  (in  the  absence  of 
desorption).  Barrier  controlled  adsorption  refers  to  a  situation  in  which  a  barrier  exists  to 
transfer  of  amphiphiles  from  the  subsurface  to  the  surface.  In  this  case  the  subsurface  and 
the  bulk  concentrations  of  amphiphiles  are  equal  and  the  kinetics  of  monolayer  formation 
can  be  characterized  by  a  well  chosen  rate  expression. 

To  understand  PC  monolayer  formation  at  the  aqueous:CCl4  interface,  we  must  first 
determine  whether  monolayer  formation  follows  diffusion  controlled  or  barrier  controlled 
kinetics.  Implicit  in  this  analysis  is  the  assumption  that  the  measured  dynamic  surface 
tension  can  be  related  to  the  time  dependent  surface  concentration  in  a  well  defined  manner. 
The  choice  of  an  appropriate  equation  of  state  relating  measured  surface  tension  to  surface 
concentration  is  addressed  a  subsequent  section. 

The  first  physically  based  model  used  to  describe  diffusion  controlled  kinetics  was 
developed  Ward  and  Tordai.45  Their  approach  drew  from  the  well  developed  theory  for 
heat  transfer  and  expressed  adsorption  (or  surface  concentration)  in  the  following  form: 

r(t)  =  2^Co^t  -  ff  c(0,  t  -  T)cWt)  ( 1 ) 

where 

T  =  surface  concentration  (mole/m2) 

c0  =  bulk  surfactant  concentration  (mole/m3) 
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D  =  Diffusion  coefficient  (m2/s) 
t  =  time  (s) 
n  =  3.1415.  .  . 

This  form  of  the  surface  concentration  does  not  possess  an  analytical  solution  but  it  can  be 
simplified  depending  on  whether  equilibration  of  the  monolayer  occurs  quickly  or  slowly 
relative  to  the  characteristic  diffusion  time  of  the  system. 48  Provided  that  the  monolayer 
remains  expanded  (low  r)>  the  surface  concentration  can  be  directly  related  to  the  measured 
surface  tension.  The  short  time  approximation  of  Eq.  1  states  that  the  surface  concentration 
(and  the  surface  tension)  scales  as  im  while  the  long  time  approximation  predicts  that 
surface  tension  should  scale  as  t'1/2:  48<51 


Short  Time  Approximation: 


(2) 


Long  Time  Approximation: 


(3) 


where 

7(t)  =  time  dependent  surface  tension  (mN/m). 

If  monolayer  formation  follows  diffusion  controlled  kinetics,  then  plotting  *y(t)  against 
either  tm  or  t'm  should  result  in  a  straight  line.  If  monolayer  formation  results  from  a 
barrier  controlled  process,  the  time  dependent  surface  concentration  (and  surface  tension) 
will  be  determined  by  an  appropriate  rate  equation. 
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Figure  4  shows  the  dynamic  surface  tension  of  the  aqueous:CCl4  interface  for  a 
1  jxMolar  solution  of  DLPC  vesicles.  The  top  two  panels  plot  Y(t)  versus  t'l/2  and  t1/2, 
respectively  (Figures  4a  and  4b).  If  monolayer  formation  were  diffusion  controlled  in  the 
long  time  limit,  then  the  tension  data  in  Figure  4a  should  exhibit  linear  behavior  closer  to 
the  origin.  Diffusion  limited  behavior  in  the  short  time  limit  should  show  7(t)  depending 
linearly  on  t1/2  in  Figure  4b.  Clearly,  7(t)  is  not  linear  in  either  the  Figure  4a  or  4b.  A  plot 
of  the  surface  tension  versus  time  appears  in  the  bottom  panel  of  Figure  4.  In  general  y(t) 
follows  a  single  exponential  decay  for  bulk  PC  concentrations  <  1  pMolar.  At  higher  ... 
concentrations  the  time  dependent  surface  tension  decays  in  a  more  complicated  fashion. 
Taken  together,  the  three  castings  of  the  dynamic  surface  tension  shown  in  Figure  4 
strongly  suggest  that  PC  monolayer  formation  at  an  aqueous:CCl4  interface  from  a  solution 
of  PC  vesicles  represents  a  barrier  controlled  rather  than  a  diffusion  controlled  process. 

Characterizing  the  time  dependent  surface  tension  for  a  barrier  controlled  process 
requires  formulating  a  rate  expression  for  the  time  dependent  surface  concentration  (n(t)) 
and  relating  n(t)  to  7(t)  by  means  of  a  2-dimensional  equation  of  state.  In  creating  a  kinetic 
model  we  assume  that  the  barrier  to  monolayer  formation  arises  from  either  extracting  PC 
monomers  from  vesicles  or  complete  vesicle  rapture  at  the  aqueous:CCl4  interface.  These 
two  mechanisms  represent  the  extraction  mechanism  and  rapture  mechanism  described 
below.  We  also  postulate  that  at  low  monolayer  concentrations  vesicles  are  not  surface 
active,  a  reasonable  assumption  given  that  the  hydrophobic  acyl  chains  of  the  PC 
monomers  lie  buried  in  the  vesicle  bilayer  interior  and  that  the  outer  surface  of  the  vesicle 
consists  of  zwitterionic  charged  headgroups.  At  higher  monolayer  concentrations 
Coulombic  attractions  between  the  headgroups  of  adsorbed  monomers  and  the  charged 
vesicle  shells  may  lead  to  higher  vesicle  concentration  at  the  interface.14  Finally,  we 
propose  that  at  sufficiently  low  PC  concentrations,  the  lopsided  partitioning  of  PC 
monomers  into  bilayers  (vs.  free  in  solution),  renders  the  effect  of  free  monomer 
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adsorption  inconsequential.  Based  on  the  partitioning  results  discussed  in  the  Experimental 
section,  we  expect  a  vanishingly  small  aqueous  phase  concentration  of  free  monomers. 


B.  Kinetic  Mechanisms 

1.  Extraction  Mechanism 

The  monomer  extraction  model  requires  that  each  PC  vesicle  deposit  some  fraction 
of  its  monomers  at  the  interface  one  monomer  at  a  time.  This  mechanism  assumes  that  the 
vesicles  themselves  remain  intact  and  that  the  number  of  vesicles  remaining  after  the 
monolayer  has  equilibrated  is  approximately  equal  to  the  number  of  vesicles  prior  to  the 
start  of  monolayer  formation.  In  order  for  this  condition  to  be  true  the  total  number  of  PC 
molecules  in  aqueous  solution  must  be  much  greater  than  the  number  of  molecules  to  form 
a  tightly  packed  monolayer  at  the  aqueous:CCl4  interface.  Otherwise,  formation  of  a 
monolayer  will  necessarily  result  in  destruction  of  vesicles.  With  this  model  the  rate  of 
monomer  adsorption  at  the  interface  can  be  described  by  simple  pseudo-zeroth  order 
kinetics: 


dn(t) 

dt 


=  kv0 


(4) 


Here,  n(t)  represents  the  time  dependent  surface  concentration,  k  equals  the  rate  constant 
for  monomer  extraction  from  the  vesicle  bilayer  and  v0  reflects  the  vesicle  concentration  in 
the  subsurface  layer  during  the  monolayer  formation  process.  Because  the  right  hand  side 
of  Eq.  4  is  a  constant  for  a  given  temperature  and  PC  concentration,  the  differential 
equation  is  readily  integrated  to  provide  a  solution  for  n(t): 

n(t)  =  kv0t  (5) 

At  sufficiently  low  concentrations,  the  monomer  extraction  mechanism  either  limits  the 
equilibrium  surface  concentration  or  begins  to  break  down  as  the  vesicles  become  too  small 
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to  remain  thermodynamically  stable,  which,  in  turn,  would  lead  to  vesicle  rupture. 

Complete  vesicle  decomposition  represents  the  second  proposed  mechanism  of  monolayer 
formation. 

2.  Rupture  Mechanism 

The  vesicle  rupture  mechanism  involves  a  PC  vesicle  approaching  the  aqueous:CCl4 
interface,  breaking  apart,  and  allowing  all  of  its  monomers  to  adsorb  to  the  interface.  The 
kinetics  of  this  model  can  be  described  by  the  following  expression: 


dn(t)  _  dv(t) 
~d T"  m  dt 


(6) 


where  n(t)  and  v(t)  are  the  time  dependent  concentrations  of  adsorbed  monomers  and 
vesicles  in  the  subsurface  layer,  respectively;  m  refers  to  the  average  number  of  monomers 
in  a  vesicle. 

In  the  low  concentration  limit  the  vesicle  concentration  depends  on  extent  of 
monolayer  formation  and  follows  first  order  kinetics: 


Here,  k  represents  the  rate  constant  for  vesicle  rupture  at  the  interface.  Conservation  of 
molecules  requires  that  the  number  of  adsorbed  monomers  equals  the  number  of  ruptured 
vesicles  multiplied  by  the  number  of  monomers  per  vesicle  (assumed  to  be  a  constant): 


n(t)  =  m(v0  -  v(t))  =  mv0(l  -  e~kt) 


(8) 
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According  to  Eq.  8,  the  time  dependent  surface  concentration  should  follow  first  order 
kinetics,  asymptotically  approaching  its  equilibrium  limit  according  to  the  rate  constant  for 
vesicle  decomposition,  k.  The  term  mv0  simply  reflects  the  initial  bulk  PC  concentration 
(which  is  equal  to  the  initial  PC  concentration  in  the  subsurface  layer). 

In  the  high  concentration  limit  the  vesicle  concentration  remains  unchanged  and  the 
solution  to  Eq.  6  for  n(t)  reduces  to  pseudo  zeroth  order  kinetics  as  in  Eq.  5. 

These  two  different  expressions  for  n(t)  arising  from  the  rupture  mechanism  .. 
demonstrate  that  the  total  bulk  phosphatidylcholine  concentration  can  play  an  important  role 
in  the  rate  of  monolayer  formation.  Experimental  isotherms  for  these  PC  systems  above 
their  respective  transition  temperatures  (Fig.  3b,  DLPC)  show  that  tightly  packed 
monolayers  consist  of  a  total  of  ~5  x  1015  monomers  (in  27  cm2)  corresponding  to 
molecular  areas  of  55  ±  9  A2/molecule.  For  25  ml  of  a  0.5  p.Molar  PC  solution,  -70%  of 
the  molecules  in  solution  would  be  required  to  form  a  tightly  packed  monolayer.  This 
percentage  does  not  take  into  account  PC  monomers  which  adsorb  to  the  walls  of  the  dish 
or  to  the  aqueousrair  interface.  Should  these  effects  prove  non-neglible,  then  higher  bulk 
concentrations  would  be  required  in  order  to  form  tightly  packed  monolayers.  In  fact  we 
do  not  observe  formation  of  tightly  packed  monolayers  from  aqueous  solutions  of  PC 
vesicles  in  their  liquid  crystalline  phase  until  bulk  concentrations  of  ~2pMolar.  At  this  bulk 
concentration,  only  -15%  of  the  PC  molecules  in  solution  need  to  adsorb  to  the 
aqueous:CCl4  interface  in  order  to  form  a  monolayer  with  the  observed  surface  coverage. 

C.  Relating  n(t)  to  y(t):  Equations  of  State 

1.  Ideal  Equation  of  State 

Relating  the  time  dependent  surface  concentration  (Eq.  5  or  Eq.  8)  to  the 
experimentally  measured  dynamic  surface  tension,  y(t),  requires  choice  of  an  appropriate 
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two  dimensional  equation  of  state.  In  the  simplest  case  the  measured  surface  tension  can  be 
related  to  the  surface  concentration  by  the  two  dimensional  ideal  gas  equation: 28 

FLA.  =  n(t)RT  (9) 

where  n  =  y0-y(i) 

A  =  interfacial  area  (m2) 

n(t)  =  surface  concentration  (mole) 

R  =  8.314  Jmole'1K'1 
T  =  temperature  (K) 

The  ideal  equation  of  state  assumes  that  the  adsorbed  molecules  are  non-interacting  and 
present  at  low  surface  concentrations. 

2.  Frumkin-based  Equation  of  State 

As  monolayers  become  more  concentrated  adsorbed  monomers  may  begin  to 
experience  intermolecular  interactions  leading  to  a  breakdown  of  Eq.  9.  An  equation  of 
state  based  on  the  Frumkin  isotherm  takes  these  forces  into  account: 28 

HA  =  n°°RT  lnfl  -  ^  (10) 


Here,  n°°  represents  the  surface  concentration  of  a  tightly  packed  monolayer. 

Substituting  the  expression  for  n(t)  from  either  kinetic  model  (Eq.  5,  Eq.  8)  into 
these  different  equations  of  state  allows  one  to  obtain  an  analytical  means  of  fitting  the 
dynamic  surface  tension  as  a  function  of  time.  Table  2  summarizes  the  results  as  Equations 
11-13.  Equations  are  grouped  (e.g.  Eq.  1  la  and  1  lb)  according  to  their  analytic 
functional  form.  Several  sets  of  conditions  predict  that  the  surface  tension  should  depend 
on  time  in  a  linear  fashion  (Eq.  1  la  and  1  lb) :  a)  the  extraction  mechanism  coupled  with 
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the  ideal  equation  of  state;  b)  the  rupture  mechanism  in  the  high  concentration  limit  coupled 
with  the  ideal  equation  of  state;  and  c)  the  rupture  mechanism  in  the  low  concentration  limit 
coupled  with  the  Frumkin-based  equation  of  state.  Assuming  the  low  concentration  limit  of 
the  rupture  mechanism  and  the  ideal  equation  of  state  predicts  that  the  surface  tension 
should  follow  an  exponential  decay  with  time  (Eq.  12).  Lastly,  when  either  the  extraction 
mechanism  or  the  rupture  mechanism  (high  concentration  limit)  is  coupled  to  the  Frumkin 
based  equation  of  state,  y(t)  assumes  a  time  dependence  which  is  not  easily  represented 
(Eq.  13a).  However,  this  expression  can  be  inverted  leading  to  a  simple  exponential  .. 
dependence  of  time  on  the  measured  surface  tension.  (Eq.  13b  and  13c). 

IV.  Results  and  Discussion 

A.  Concentration  Dependence  of  Monolayer  Formation 

1.  Low  Concentration  Limit 

Experiments  carried  out  above  the  bilayer  gel-liquid  crystalline  transition 
temperature  (Tc)  show  the  rate  and  extent  of  monolayer  formation  to  depend  quite 
sensitively  on  the  bulk  PC  concentration.  The  DLPC  isotherm  in  Fig.  3b  demonstrates 
how  the  equilibrium  properties  of  the  monolayer  change  with  bulk  PC  concentration.  At 
bulk  concentrations  of  >  2  fiMolar,  the  terminal  surface  pressure  reaches  an  asymptotic 
limit  of  ~42  mN/m.  Further  increase  in  bulk  concentration  does  not  lead  to  any  additional 
increase  in  surface  pressure,  indicating  an  upper  limit  to  the  extent  of  monomer  packing  at 
the  aqueous:CCl4  interface.  As  mentioned  in  the  experimental  section,  analysis  of  the 
DLPC  isotherm  data  by  means  of  the  Gibbs  equation  calculates  that  these  tightly  packed 
monolayers  have  concentrations  of  1.8  x  1014  molecules/cm2  or,  equivalently,  molecular 
areas  of  55  A2/molecule.30*31 

Although  the  equilibrium  or  terminal  surface  pressure  does  not  change  for  bulk  PC 
concentrations  above  2  (iMolar,  the  rate  at  which  the  monolayers  form  still  depends  upon 
bulk  concentration.  The  upper  panel  of  Figure  5  shows  the  time  dependent  surface  tension 
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during  monolayer  formation  from  solutions  of  different  bulk  DLPC  concentrations.  At 
bulk  concentrations  of  ~1  jxMolar  and  below,  the  data  decay  in  an  exponential  fashion 
slowly  approaching  an  asymptotic  limit  determined  by  bulk  concentration  and  experimental 
temperature.  The  lower  panel  of  Figure  5  shows  a  expanded  view  of  the  ljiMolar  system. 
Indistinguishable  from  the  data  at  all  but  the  earliest  times  (<  10  min)  is  the  fit  to  an 
exponential  decay.  Except  for  the  earliest  times  during  monolayer  formation,  the 
exponential  decay  reproduces  the  experimental  data  to  within  0. 1  mN/m,  well  within  the 
limits  of  experimental  uncertainty.  This  behavior  of  the  surface  tension  is  consistent  with 
Equation  12  which  couples  the  vesicle  rupture  mechanism  to  the  ideal  equation  of  state. 

Fitting  the  data  to  a  single  exponential  decay  requires  that  the  fit  coefficients  assume 
physically  meaningful  definitions.  Equation  12  has  the  general  form: 

y(t)  =  ao  +  (14) 

From  Eq.  12  we  see  that  ^  corresponds  to  (y0  -  RTmvc/A)  and  a,  equals  RTmv0/A.  The 
term  in  the  exponent,  is  simply  the  rate  constant  for  vesicle  rupture,  k.  In  physical 
terms,  ao  should  equal  the  equilibrium  interfacial  tension  for  a  system  with  a  specified 
DLPC  bulk  concentration,  a,  should  scale  linearly  with  bulk  concentration  (according  to  the 
mv0  term)  and  %  is  not  expected  to  show  any  concentration  dependence. 

Figure  6  plots  the  fit  parameters,  ao,  a,  and  a^  against  bulk  DLPC  concentration. 

At  concentrations  below  2  jiMolar,  the  fit  parameters  behave  in  the  manner  predicted  by 
Eq.  12.  Fitting  the  2  pMolar  system  leads  to  a  physically  unreasonable  ao  value  of  -5 
mN/m  while  the  ao  parameter  for  higher  concentrations  varies  without  a  systematic  trend. 
The  a,  coefficient  depends  linearly  on  bulk  DLPC  concentration  up  to  a  bulk  concentration 
of  1 .2  pMolar.  At  concentrations  of  2  pMolar  and  higher  the  a,  coefficient  levels  out  at  a 
value  of  ~30  mN/m.  We  note  that  at  the  lower  concentrations,  the  a,  coefficients  from 
calculated  fits  correlate  reasonably  well  with  the  value  one  would  calculate  using  the 
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analytical  expression  in  Eq.  12  (a,  =  RTmv/A).  For  example,  at  a  bulk  concentration  of 
1  (iMolar,  the  analytical  form  of  a,  calculates  a  value  of  20  mN/m  while  the  fit  yields  a 
value  of  15  mN/m.  Given  that  the  model  does  not  take  into  account  loss  of  PC  monomers 
to  the  air:aqueous  interface  or  to  the  walls  of  the  dish,  the  “effective”  PC  concentration 
(mv0)  may  very  well  be  lower  than  the  reported  concentration.  This  overestimation  of  the 
bulk  PC  concentration  by  the  analytical  expression  would  lead  to  a  calculated  a,  value 
which  was  higher  than  the  value  from  a  fit  to  Eq.  14  which  samples  the  experimental  or 
effective  PC  concentration.  Finally,  the  plot  of  versus  concentration  shows  the  rate  .„ 
constant  for  monomer  adsorption  to  be  constant  at  sufficiently  low  concentrations 
(<  2pMolar).  At  higher  concentrations,  the  rate  constant  appears  to  climb  but  this 
observation  should  be  tempered  by  the  fact  that  at  these  higher  concentrations,  Eq.  12  no 
longer  captures  the  functional  form  of  the  dynamic  surface  tension. 

2.  High  Concentration  Limit 

At  bulk  PC  concentrations  >  2  pMolar  the  DLPC  isotherm  (Figure  3b)  indicates  that 
the  resulting  monolayer  approaches  its  tight  packed  limit.  Consequently,  the  two 
dimensional  ideal  equation  of  state  (Eq.  9)  should  no  longer  accurately  describe  the 
adsorbed  monolayer.  More  appropriate  is  an  equation  of  state  such  as  Eq.  10  which  takes 
into  account  intermolecular  interactions  between  monomers  at  the  interface.  Because 
experimental  conditions  already  specify  a  high  concentration  of  PC  vesicles  in  bulk 
solution,  the  expression  for  the  time  dependent  surface  concentration,  n(t),  does  not  depend 
on  the  choice  of  mechanism.  Both  the  extraction  and  the  rupture  mechanisms  express  n(t) 
as  a  linear  function  of  time  (Eq.  5).  Substituting  into  Equation  10  for  n(t)  leads  to  an 
expression  of  y(t)  in  terms  of  time  which  is  not  easily  represented  (Eq.  13a.,  Table  2). 
However,  Eq.  13a  may  be  inverted  to  recast  time  as  a  function  of  interfacial  tension,  y. 
Carrying  out  this  operation  leads  to  Eq.  13b/c  in  Table  2. 
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Figure  7  shows  the  results  of  trying  to  fit  the  high  concentration  surface  tension 
data  with  Eq.  12  and  Eq.  13b/c.  By  inspection  we  can  eliminate  Eq.  1 1  which  predicts 
linear  dependence  of  y(t)  on  time.  In  the  top  panel  of  Fig.  7,  surface  tension  is  plotted  vs. 
time  and  the  data  are  fit  to  a  single  exponential  decay  according  to  Eq  12.  Although  the 
differences  between  calculated  and  experimental  surface  tensions  never  grow  larger  than 
1.5  mN/m,  close  analysis  reveals  that  the  fit  deviates  from  the  data  in  a  systematic  fashion 
and  never  accurately  describes  the  time  dependent  behavior  of  the  surface  tension.  In  the 
bottom  panel  the  experimental  time  is  plotted  versus  the  measured  surface  tension  and  tlie 
data  is  fit  according  to  Eq.  13b/c.  Between  surface  tensions  of  43  mN/m  (t  =  0)  and  6 
mN/m  (t  =  120  min)  the  fit  is  virtually  indistinguishable  from  the  data.  We  take  the 
remarkable  correlation  between  the  experimental  data  and  the  calculated  fit  as  an  indication 
that  Equation  13b/c  does,  in  fact,  accurately  capture  the  important  aspects  of  the  monolayer 
formation  process  in  the  high  aqueous  PC  concentration  limit.  Furthermore,  if  one 
assumes  the  vesicle  rupture  mechanism,  then  the  rate  constant,  k  (=n“/a0mv0)  agrees  with 
vesicle  rupture  rate  constants  taken  from  the  exponential  term  in  the  low  concentration 
experiments.  This  high  concentration  rate  constant  appears  as  the  open  triangle  in  Fig.  5c. 

After  the  first  two  hours  of  the  8  pMolar  experiment,  the  surface  tension  slowly 
decays  from  6  mN/m  to  4  mN/m  in  a  manner  inconsistent  with  the  analysis  between 
0  min  <  t  <  120  min.  This  behavior  suggests  a  crossover  in  the  mechanism  of 
monolayer  formation,  possibly  from  the  proposed  rupture  mechanism  to  more  of  an 
extraction  type  mechanism.  Intuitively,  such  a  crossover  seems  reasonable  given  that 
during  the  latter  stages  of  monolayer  formation  at  high  PC  concentrations,  vesicles  no 
longer  approach  a  neat  aqueous:CCl4  interface.  Rather,  vesicles  see  a  moderately  packed 
sheet  of  PC  headgroups,  and  the  limiting  constraint  on  the  monolayer  formation  kinetics 
may  simply  become  creating  space  to  accommodate  the  remaining  monomers. 

Although  Figure  7  supports  using  the  Frumkin  based  equation  of  state  (Eq.  10)  for 
describing  monolayer  formation  kinetics  in  the  high  concentration  limit,  the  analysis  can  not 
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distinguish  which  mechanism  leads  to  monolayer  formation.  At  high  concentrations  the 
extraction  and  the  rupture  mechanisms  predict  a  linear  time  dependence  for  n(t)  (Eq.  5)  and, 
consequently,  similar  expressions  for  y(t)  (Eq.  11,  Eq.  13b  and  13c).  The  only  difference 
between  Eq.  13b  and  Eq.  13c  lies  in  a  factor  of  m,  the  number  of  monomers  per  vesicle 
which  we  assume  to  be  an  ensemble  averaged  constant.  However,  we  believe  that  the 
rupture  mechanism  represents  the  primary  source  by  which  monolayers  form  in  the  high 
concentration  limit  based  on  data  from  the  low  concentration  experiments.  At  low 
concentrations  only  the  rupture  mechanism  coupled  with  the  ideal  equation  of  state  lead}?  to 
an  analytic  expression  (Eq.  12)  which  accurately  fits  the  experimental  dynamic  surface 
tension  data.  A  priori  we  see  no  reason  to  assume  that  the  mechanism  during  the  initial 
stages  of  monolayer  formation  should  change  with  increasing  PC  concentration. 

Agreement  of  the  one  adjustable  parameter  in  Eqs.  12  and  13c,  the  rate  constant  for  vesicle 
rupture  k,  strengthens  this  claim.  We  therefore  propose  that  monolayer  formation  initially 
proceeds  by  means  of  vesicle  rupture.  Only  in  the  latter  stages  of  monolayer  formation 
from  solutions  having  high  PC  concentrations  do  we  expect  the  rupture  mechanism  to  give 
way  to  an  extraction  type  of  mechanism. 

B.  Dependence  of  Monolayer  Formation  on  PC  bilayer  Phase 

The  previous  section  provided  evidence  supporting  the  rupture  mechanism  as  the 
means  by  which  PC  monolayers  form  at  an  aqueous:CCl4  interface  from  a  solution  of 
aqueous  phase,  liquid  crystalline  vesicles.  Dynamic  surface  tension  experiments  in  both 
the  low  and  high  concentration  limits  were  consistent  with  predictions  of  the  rupture 
mechanism  coupled  with  the  appropriate  equations  of  state.  We  now  show  how  both  the 
rate  and  extent  of  monolayer  formation  depends  sensitively  on  lipid  bilayer  phase.  Figure  8 
contains  time  dependent  surface  tension  data  for  lpMolar  solutions  of  DLPC  (n  =  12), 
DMPC  (n  =  14),  DPPC  (n  =  16)  and  DSPC  (n  =  18)  at  22°  C.  Shown  are  the  first  10  hrs 
of  each  experiment  as  well  as  the  asymptotic  surface  tension  limits  for  each  system.  The 
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time  required  to  reach  these  limits  varies  from  ~24  hrs  (DLPC)  to  >  72  hrs  (DSPC  and 
DPPC  systems).  Three  days  marked  the  limit  of  our  ability  to  collect  reliable  data  from  a 
single  experiment.  After  three  days  the  age  of  the  interface  became  a  concern  and 
evaporative  losses  began  to  play  a  role  in  the  surface  tension  measurements. 

The  data  in  Fig.  8  clearly  show  that  at  room  temperature  the  extent  of  monolayer 
formation  depends  inversely  on  the  hydrocarbon  chain  length  of  the  PC  species.  DLPC 
forms  the  most  concentrated  monolayer,  while  DPPC  and  DSPC  form  the  most  expanded 
monolayers.  Transition  temperatures  in  Table  1  show  that  for  these  ambient  experiments 
(T  =  22°  C),  the  DLPC  solution  contains  vesicles  in  their  more  fluid,  liquid  crystalline 
phase.  DMPC  vesicles  exist  in  a  liquid  crystalline/gel  coexistence  region  while  DPPC  and 
DSPC  vesicles  are  frozen  with  their  lipid  bilayers  in  a  rigid  gel  phase.  Despite  these 
differences  the  systems  do  share  similarities.  Most  notably,  surface  tension  experiments 
show  evidence  that  monolayer  formation  remains  barrier  controlled  regardless  of  vesicle 
bilayer  phase,  and  kinetics  in  the  low  concentration  limit  always  appear  to  be  first  order  in 
bulk  PC  concentration. 

Both  the  DSPC  isotherm  (Fig  3b.)  and  the  time  dependent  surface  tension  data 
suggest  that  PC  vesicle  solutions  below  Tc  can  only  form  expanded  monolayers  and  the 
monolayers  themselves  form  very  slowly.  In  earlier  work  we  demonstrated  how  the  extent 
of  monolayer  formation  changes  discontinuously  at  the  gel/liquid  crystalline  transition 
temperature.31  At  temperatures  above  their  respective  transition  temperatures,  PC  solutions 
with  concentrations  £  2  jiMolar  irreversibly  form  monolayers  with  surface  coverages  of 
55  A2/molecule.  Below  Tc  PC  solutions  form  expanded  monolayers  regardless  of  bulk 
concentration  provided  that  the  monolayers  are  allowed  to  form  without  external 
perturbation.  Recent  work  in  this  lab  has  shown  that  physically  disturbing  the  interface 
following  partial  monolayer  formation  can  produce  compact  monolayers  even  at 
temperatures  below  Tc.52  Under  such  circumstances  monolayers  of  longer  chain  species 
exhibit  characteristics  consistent  with  very  tightly  packed  monomers. 
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C.  Temperature  Dependence  of  Monolayer  Formation 

The  dynamic  surface  tension  experiments  suggest  fundamental  differences  in  the 
monolayer  formation  process  depending  on  whether  the  experimental  temperature  lies 
above  or  below  Tc.  Although  low  concentration  data  support  the  rupture  mechanism 
(Eq.  12)  regardless  of  bilayer  phase,  rupture  of  a  liquid  crystalline  vesicle  may  proceed  in 
a  very  different  manner  than  rupture  of  a  gel  phase  vesicle.  In  order  to  examine  the  effects 
of  the  lipid  bilayer  phase  on  monolayer  formation,  we  measured  the  dynamic  surface 
tension  of  the  four  PCs  in  the  low  concentration  limit  at  different  temperatures. 

Figure  9  shows  representative  dynamic  surface  tension  data  for  0.5  pMolar  DLPC 
solutions  as  a  function  of  temperature.  Included  on  the  graph  are  the  exponential  decays 
which  result  from  fitting  the  data  to  Eq.  12.  All  four  experiments  exhibit  exponential  decay 
behavior,  and  from  the  fit  parameters  emerges  the  temperature  dependence  of  the 
monolayer  formation  rate.  Plotting  the  natural  log  of  the  rate  constant  (a^  against  inverse 
temperature  generates  an  Arrhenius  plot  from  which  an  activation  energy  for  monolayer 
formation  may  be  calculated. 

Arrhenius  plots  for  all  four  phosphatidylcholines  appear  in  Figure  10  and  the  results 
are  summarized  in  Table  3.  The  liquid  crystalline  DLPC  data  (Fig.  10a)  lead  to  a 
monolayer  formation  activation  energy  (Et)  of  39  ±  7  kJ/mol.  Activation  energies  for  both 
gel  phase  DPPC  (Fig.  10c)  and  gel  phase  DSPC  (Fig.  lOd)  are  considerably  lower,  17  ±  6 
kJ/mol  and  20  ±11  kJ/mol,  respectively.  The  CI4  species,  DMPC,  affords  a  unique 
opportunity  to  explicitly  examine  how  the  phase  of  the  lipid  bilayer  affects  the  monolayer 
formation  process.  With  a  Tc  of  23°  C  the  temperature  dependence  of  monolayer  formation 
kinetics  can  be  examined  for  solutions  of  vesicles  both  in  the  rigid  gel  phase  (T  <  23°  C) 
and  in  the  more  fluid,  liquid  crystalline  phase  (T  >  23°  C).  The  DMPC  kinetic  data  (Fig. 
10b)  show  a  break  in  the  Arrhenius  plot  at  approximately  the  DMPC  transition  temperature. 
Monolayer  formation  from  aqueous  solutions  of  liquid  crystalline  vesicles  has  an  E0  of  41  ± 
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8  kJ/mol,  effectively  identical  to  that  of  DLPC.  However,  for  solutions  of  vesicles  in  their 
gel  state,  Ea  drops  by  a  factor  of  five  to  8  ±  5  kJ/mol.  Generally  speaking,  Ea  appears  to 
be  larger  for  solutions  of  liquid  crystalline  vesicles  than  for  solutions  of  gel  phase  vesicles. 

While  studying  rates  of  phosphatidylcholine  monomer  exchange  between 
unilamellar  vesicles,  McLean  and  Phillips  observed  similar,  lipid  bilayer  phase  dependent 
behavior.22  Specifically,  exchange  kinetics  of  DMPC  monomers  between  DMPC  vesicles 
in  their  liquid  crystalline  phase  exhibited  a  much  stronger  temperature  dependence  than  did 
monomer  exchange  between  gel  phase  vesicles.  The  calculated  activation  energies  for  :j 
monomer  exchange  were  70  kJ/mol  for  exchange  between  liquid  crystalline  vesicles  and 
47  kJ/mol  for  gel  phase  vesicles  22  Although  both  of  these  values  lie  considerably  higher 
than  the  calculated  activation  energies  for  monolayer  formation  (41  kJ/mol  liquid 
crystalline,  8  kJ/mol  gel)  the  trend  of  a  larger  barrier  for  liquid  crystalline  systems  is 
consistent  with  the  results  of  the  temperature  dependent  dynamic  surface  tension 
experiments.  Quantitative  differences  in  the  energetics  of  the  two  processes  -  monomer 
exchange  and  monolayer  formation  -  may  arise  due  to  the  different  mechanisms  by  which 
the  two  phenomena  proceed.  Monomer  exchange  between  vesicles  necessarily  requires  an 
extraction/insertion  type  mechanism20  while  the  dynamic  surface  tension  experiments 
discussed  in  the  previous  sections  suggest  that  monolayer  formation  proceeds  via  vesicle 
rupture.  If  monomer  insertion  into  a  host  bilayer  has  a  large  energetic  barrier  relative  to 
adsorption  at  the  aqueous:CCl4  interface,  then  one  would  expect  monomer  exchange  to 
have  higher  activation  energies  than  monolayer  formation. 

Despite  differences  between  the  mechanisms  of  monomer  exchange  and  monolayer 
formation,  we  expect  the  concepts  used  to  explain  the  phase  dependent  activation  energies 
to  transfer  between  the  two  phenomena.  Based  on  earlier  work  carried  out  on  micellar 
dynamics,53  McLean  and  Phillips  postulated  that  the  smaller  barrier  for  exchange  between 
gel  phase  vesicles  arose  from  the  high  degree  of  order  amongst  the  acyl  chains  within  the 
bilayer.22  A  phospholipid  monomer  contained  within  a  gel  phase  bilayer  already  has  its 
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chains  aligned  in  all -trans  conformations  and  can  easily  slide  out  of  the  bilayer  into  solution 
given  a  strong  enough  chemical  driving  force.  4,5,54  jn  contrast  acyl  chains  of  a 
phospholipid  monomer  in  a  liquid  crystalline  bilayer  are  considerably  more  disordered  and 
require  straightening  before  the  monomer  can  separate  from  the  vesicle  wall.54-55  (See,  for 
example,  Fig.  2b.)  This  additional  energy  needed  to  straighten  the  chains  of  the  departing 
monomer  leads  to  the  higher  activation  energy  for  monomer  exchange  between  liquid 
crystalline  vesicles  53 

By  analogy,  we  speculate  that  gel  phase  vesicles  face  a  smaller  barrier  to  vesicle, 
rupture  at  the  aqueous:CCl4  interface  due  to  the  well  ordered,  all -trans  hydrocarbon  chains 
which  make  up  the  bilayer  interior.  The  rigid,  gel  phase  vesicles  will  be  more  susceptible 
to  vesicle  fracture  than  the  more  fluid  liquid  crystalline  vesicles  whose  bilayer  interiors 
possess  a  much  higher  degree  of  disorder.  In  liquid  crystalline  vesicles  acyl  chain 
entanglement  and  flexibility  should  inhibit  vesicle  rupture  at  the  interface  which  results  in  a 
higher  activation  energy  for  monolayer  formation.  We  note  that  if  the  nature  of  vesicle 
rupture  changes  as  vesicles  pass  from  the  gel  state  to  the  liquid  crystalline  state,  we  might 
expect  to  see  this  effect  appear  in  the  rate  of  monolayer  formation  at  Tc.  The  DMPC  data  in 
Figure  9b  provide  evidence  of  such  a  change  with  the  rate  constant  for  monolayer 
formation  just  below  Tc  being  slightly  -  but  reproducibly  -  larger  than  the  rate  constant  just 
above  Tc.  At  Tc  monolayer  formation  ffom  gel  phase  vesicles  should  proceed  more  quickly 
than  ffom  liquid  crystalline  vesicles  due  to  the  differences  in  the  monolayer  formation 
activation  energies. 

The  disparity  Ea  for  gel  vs.  liquid  crystalline  phase  vesicles  prompts  the  question  of 
why  gel  phase  vesicles  with  a  smaller  barrier  to  monolayer  formation  only  form  expanded 
monolayers  regardless  of  bulk  PC  concentration,  while  liquid  crystalline  vesicles  form 
tightly  packed  monolayers  at  bulk  PC  concentrations  above  2  (iMolar.  Resolving  this 
paradigm  requires  treating  monolayer  formation  as  a  thermodynamically  driven  process 
where  energetics  of  initial  states  (monomers  in  vesicles)  and  final  states  (monomers  in 
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monolayers)  must  be  considered.  Figure  1 1  depicts  a  schematic  diagram  illustrating  the 
monolayer  formation  process.  The  vesicle  rupture  reaction  coordinate  appears  on  the 
abscissa  while  the  ordinate  represents  the  internal  energy  of  the  reactants  and  the  products. 
When  a  system  of  vesicles  lies  just  above  and  just  below  Tc,  liquid  crystalline  vesicles  are 
energetically  separated  from  gel  phase  vesicles  by  AHfus  for  the  gel-liquid  crystalline  phase 
transition,  a  quantity  measured  in  differential  scanning  calorimetry  (DSC) 
experiments.56’57  DSC  data  for  the  four  phosphatidylcholines  studied  show  AH^  to  be 
17  kJ/mol  for  DLPC,  26  kJ/mol  for  DMPC,  40  kJ/mol  for  DPPC,  and  45  kJ/mol  for _ 
DSPC.56’57 

Appearing  between  the  initial  and  final  states  of  the  vesicle/interface  system  are  the 
barriers  to  gel  and  liquid  crystalline  vesicle  rupture.  The  Arrhenius  data  in  Figure  10  and 
Table  3  show  the  barrier  to  be  2-5  times  larger  for  solutions  of  liquid  crystalline  vesicles 
than  for  solutions  of  vesicles  in  their  gel  phase.  In  order  for  a  monolayer  to  form,  this 
barrier  must  be  overcome.  This  picture  may  oversimplify  the  actual  monolayer  formation 
process  in  the  limit  when  tightly  packed  monolayers  form  from  solutions  of  liquid 
crystalline  vesicles.  Such  circumstances  may  require  that  a  second  barrier  appear  in  order 
to  represent  adsorption  of  monomers  after  the  monolayer  has  mostly  formed.  The  high 
concentration  data  in  Figure  7b  already  hint  that  the  mechanism  changes  during  the  latter 
stages  of  monolayer  formation  (at  t  >  120  min).  However,  in  the  low  concentration  limit, 
we  assume  that  if  a  vesicle  ruptures  at  the  interface,  its  monomers  adsorb  without  having  to 
overcome  any  further  kinetic  or  thermodynamic  barriers. 

On  the  right  hand  side  of  Figure  1 1  are  different  stages  of  PC  monolayer  formation. 
As  the  monolayer  becomes  more  concentrated  (as  more  monomers  adsorb  to  the  interface) 
the  internal  energy  of  the  monolayer  rises.  We  stress  that  this  picture  does  not  attempt  to 
quantitatively  relate  the  internal  energies  of  the  gel  and  liquid  crystalline  vesicles  to  the 
internal  energies  of  the  different  monolayer  states.  Nevertheless,  we  may  draw  several 
qualitative  conclusions  from  Figure  11:  when  faced  with  a  neat  aqueous:CCl4  interface, 
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solutions  of  PC  vesicles  experience  a  driving  force  to  rupture  and  form  monolayers 
regardless  of  lipid  bilayer  phase.  Initially,  vesicle  rupture  probably  represents  an 
entropically  driven  process.  The  system  will  become  increasingly  disordered  as  the  first 
vesicles  disintegrate  depositing  their  monomers  at  the  interface.  Enthalpic  contributions  to 
vesicle  rupture  and  monomer  adsorption  probably  cancel  out  given  that  heats  of  solvation 
for  alkyl  chains  in  hydrocarbon  liquids  and  CC14  are  similar.58 

As  the  monolayer  forms,  its  internal  energy  rises  just  as  in  Langmuir  trough 
compression  experiments.  (In  Langmuir  trough  experiments,36  the  number  of  adsorbed 
monomers  remains  constant  while  the  monolayer  area  changes;  in  these  dynamic  surface 
tension  experiments,  the  area  remains  unchanged  while  the  number  of  adsorbed  monomers 
increases.)  Eventually,  the  internal  energy  of  the  adsorbed  monolayer  equals  that  of  the  gel 
phase  vesicle  and  the  driving  force  responsible  for  gel  phase  vesicle  rupture  abates. 
Vesicles  in  their  liquid  crystalline  phase,  however,  lie  higher  in  energy  thus  they  will 
continue  to  rupture  forming  an  increasingly  concentrated  monolayer. 

V.  Conclusions 

This  paper  has  presented  dynamic  surface  tension  experiments  which  examine  the 
formation  of  phosphatidylcholine  monolayers  at  an  aqueous:carbon  tetrachloride  interface 
from  solutions  of  PC  vesicles  by  means  of  the  Wilhelmy  plate  method.  Results  show  that 
monolayer  formation  represents  an  activated  rather  than  a  diffusion  controlled  process.  A 
number  of  different  variables  control  the  rate  of  formation  including  aqueous 
phosphatidylcholine  concentration,  temperature  and  vesicle  bilayer  phase. 

The  rate  at  which  these  monolayers  form  can  be  quantitatively  described  by  means 
of  a  simple  first  order  kinetic  model  which  assumes  that  vesicles  rupture  at  a  planar 
boundary  between  the  two  immiscible  liquids.  By  this  model  every  rupture  event  allows  a 
vesicle’s  monomers  to  adsorb  to  the  interface.  In  order  to  relate  the  time  dependent  surface 
concentration,  n(t),  to  the  experimentally  measured  time  dependent  surface  tension  Cy(t)), 
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the  two  quantities  must  be  related  by  an  appropriate  two  dimensional  equation  of  state.  We 
demonstrated  that  the  two  dimensional  ideal  gas  law  describes  the  expanded  monolayer 
experimental  data  with  remarkable  accuracy.  As  the  monolayer  grows  more  concentrated, 
the  ideal  gas  law  begins  to  break  down  requiring  the  use  of  a  Frumkin  based  equation  of 
state.  These  results  appeal  to  intuition:  at  low  surface  concentrations  we  would  expect  the 
adsorbed  monomers  to  be  non-interacting,  while  at  higher  surface  concentrations 
intermolecular  interactions  between  the  adsorbed  monomers  should  begin  to  affect 
monolayer  characteristics. 

Both  dynamic  surface  tension  measurements  (Fig.  7)  and  equilibrium  isotherms 
(Fig.  3b)  show  that  the  vesicle  bilayer  phase  plays  an  important  role  in  determining  the  rate 
and  extent  of  monolayer  formation.  For  experiments  carried  out  above  the  vesicle  bilayer 
gel-liquid  crystalline  transition  temperature,  Tc,  tightly  packed  monolayers  form  at  bulk 
concentrations  greater  than  2  pMolar  and  the  equilibration  time  is  comparatively  short . 
Experiments  carried  out  at  temperatures  below  Tc  lead  to  expanded  monolayers  regardless 
of  bulk  PC  concentration  and  experiments  require  days  to  approach  an  asymptotic  limit. 
Nevertheless,  monolayer  formation  remains  barrier  controlled  regardless  of  vesicle  bilayer 
phase  and  the  rate  of  monolayer  formation  still  behaves  in  a  manner  consistent  with  the 
rupture  mechanism.  These  observations  suggest  that  vesicle  rupture  happens  differently 
for  gel  and  liquid  crystalline  vesicles. 

In  order  to  examine  further  the  effects  of  vesicle  bilayer  phase  on 
phosphatidylcholine  monolayer  formation,  we  carried  out  dynamic  surface  tension 
experiments  at  different  temperatures  for  all  four  PCs.  Monolayer  formation  from 
solutions  of  liquid  crystalline  vesicles  had  a  relatively  large  activation  energy  of  ~40  kJ/mol 
while  monolayer  formation  from  gel  state  vesicles  was  considerably  lower  (8-20  kJ/mol). 
Solutions  of  DMPC  vesicles  provided  the  most  striking  example  of  this  trend:  below  Tc 
formation  of  a  DMPC  monolayer  encounters  an  energetic  barrier  of  8  kJ/mol  while  above 
Tc  the  activation  energy  climbs  by  a  factor  of  five  to  41  kJ/mol.  Other  studies  examining 
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intervesicle  monomer  exchange  have  encountered  a  similar  phenomenon,  namely  monomer 
exchange  between  gel  phase  vesicles  had  a  much  smaller  Ea  than  monomer  exchange 
between  liquid  crystalline  vesicles.  Both  the  monolayer  formation  and  the  monomer 
exchange  results  can  be  understood  in  terms  of  the  ordering  amongst  hydrocarbon  chains 
comprising  the  vesicle  bilayer.  Gel  phase  vesicles  have  well  ordered,  rigid  bilayer  interiors 
which  may  be  more  brittle  and  prone  to  rupture  in  the  energetically  anisotropic 
environments  presented  by  interfaces.  Bilayers  in  liquid  crystalline  vesicles  are  more 
disordered  and  flexible,  resulting  in  a  larger  barrier  to  vesicle  rupture.  The  propensity  of 
liquid  crystalline  vesicles  to  form  tightly  packed  monolayers  may  reflect  the  energetics  of 
the  initial  and  final  states  (vesicles  and  monolayers)  in  the  monolayer  formation  process. 

This  study  marks  the  first  attempt  to  quantitatively  describe  the  formation  of 
phosphatidylcholine  monolayers  at  a  planar  liquiddiquid  interface  from  solutions  of  PC 
vesicles.  While  the  models  in  this  paper  provide  a  consistent,  satisfactory  explanation  of 
the  data,  several  extensions  of  this  work  immediately  suggest  themselves.  To  test  the 
generality  of  the  methods  described  above,  similar  analyses  can  be  performed 
incorporating  different  two  dimensional  equations  of  state  or  new  mechanisms  of 
monolayer  formation.  Experimental  extensions  involve  tailoring  the  composition  of  the 
vesicle  bilayers  to  affect  the  rate  of  monolayer  formation  as  well  as  investigating  the 
structural  integrity  of  vesicles  after  incorporating  biological  macromolecules  into  the 
bilayer. 
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Table  1.  Properties  of  Phosphatidylcholines 


Molecule 

n 

Tc56 

DLPC 

12 

-rc 

DMPC 

14 

23°  C 

DPPC 

16 

U 

0 

DSPC 

18 

54°  C 
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Table  2.  Expressions  for  7(t)  from  extraction  and  rupture  mechanisms  coupled  with  ideal 
and  Frumkin  equations  of  state. 


Equation 

of 

State 


Ideal 


Frumkin 


Mechanism 

Extraction  Rupture 


Y(t)  =  Yo 


RTkvc 

A 


(Equation  11a) 


low  concentration 

(Equation  12) 


high  concentration 


Y(0  =  Yo 


RTkmv0 

A 


(Equation  lib) 


Y(0  =  Yo 


RTn° 


A 

(Equation  13a) 


low  concentration 

.  N  RTkmv0  , 
Y(t)  =  Yo - a-2-1 

(Equation  lib) 


nc 


kv 


(  -Yo  Y(t)  ^ 
1  -  eRTn°°  eRTn01 


o 


high  concentration 


t  = 


n° 


kmv 


o 


(  -Yo  Y(0 
1  —  eRTn°°  eRTn“ 


(Equation  13b) 


(Equation  13c) 
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Table  3:  Activation  energies  for  monolayer  formation 


Molecule 

Bilayer  phase 

Ea  (kJ/mol) 

DLPC 

liquid  crystalline 

39  ±7 

DMPC 

liquid  crystalline 

41  ±8 

gel 

8  ±  5 

DPPC 

gel 

17  +  6 

DSPC 

gel 

20±  11 

34 


Figure  Captions 

Figure  1.  Illustration  of  PC  monolayer  formation  at  an  aqueous:carbon  tetrachloride 
interface  from  a  solution  of  PC  vesicles.  Based  on  dynamic  light  scattering 
measurements,  we  believe  the  vesicles  used  in  these  studies  are  multilamellar 
rather  than  unilamellar  in  nature. 

Figure  2.  a)  Phosphatidylcholine  structure  b)  Schematic  picture  illustrating  the  gel-liquid 
crystalline  transition  in  vesicle  bilayers.  Table  1  contains  the  transition 
temperatures,  Tc,  for  the  phosphatidylcholines  studied  in  this  work.  ;J 

Figure  3.  a)  Schematic  representation  of  the  Wilhelmy  plate  method  for  measuring  the 
tension  of  the  aqueous:  carbon  tetrachloride  interface.  In  this  picture,  the 
aqueous  phase  wets  the  platinum  plate  exerting  an  upward  force  which  equals 
the  product  of  the  tension  and  the  plate  perimter.  b)  Isotherms  for  DLPC  and 
DSPC  under  ambient  conditions  (22°  C).  Solutions  of  DLPC  vesicles  form 
tightly  packed  monolayers  (1.8  x  1014  molecules/cm2  or  55  A2/molecule)  at 
concentrations  above  ~2  jiMolar.  Solutions  of  DSPC  vesicles  form  expanded 
monolayers  (as  evidenced  by  lower  surface  pressures)  regardless  of  bulk 
concentration. 

Figure  4.  a)  Dynamic  surface  tension  versus  time  (min)‘l/2.  If  monolayer  formation 

followed  diffusion  controlled  kinetics  in  the  long  time  limit,  the  data  would  be 
linear,  especially  near  the  origin.  The  data  begins  at  t  =  90  min  (far  right)  as 
marked  in  Fig.  4c.  b)  Dynamic  surface  tension  versus  time  (min)1/2.  Curvature 
in  the  data  indicate  that  monolayer  formation  does  not  follow  diffusion  controlled 
kinetics  in  the  short  time  limit,  c)  Dynamic  surface  tension  (solid  line)  versus 
time  (min)  and  the  fit  of  the  data  to  a  single  exponential  decay  (dashed  line).  The 
arrow  marks  the  starting  point  of  the  data  shown  in  Fig.  4a.  Note  that  the  data 
extends  out  to  ten  hours. 
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Figure  5.  a)  Dynamic  surface  tension  data  from  DLPC  solutions  with  varying  bulk  DLPC 
concentrations,  b)  The  first  six  hours  of  a  DLPC  dynamic  surface  tension 
experiement  in  the  low  concentration  limit  (same  data  as  in  Figure  4).  Data  were 
fit  to  a  single  exponential  decay  (dotted  line)  as  predicted  by  Equation  12  in 
Table  2. 

Figure  6.  Fit  parameters  a0,  a,,  and  from  fitting  DLPC  dynamic  surface  tension  data 
with  Equation  12.  Fit  parameters  agree  with  analytically  calculated  values 
through  bulk  DLPC  concentrations  of  ~1  jxMolar.  Also  plotted  on  the  a2  graph 
as  an  open  triangle  is  the  rate  constant,  k,  calculated  from  fitting  the  8  pMolar 
data  to  Equation  13b.  Agreement  between  the  high  concentration  rate  constant 
and  rate  constants  calculated  via  Equation  13b  and  the  low  concentration  rate 
constants  calculated  via  Equation  12  support  the  idea  that  monolayer  formation 
proceeds  via  the  rupture  mechanism. 

Figure  7.  a)  Dynamic  surface  tension  data  for  8  (iMolar  DLPC  plotted  versus  time 

according  to  Equation  12.  Although  the  exponential  decay  (dashed  line)  captures 
the  general  behavior  of  the  data,  it  systematically  misses  the  data  at  the  points  of 
steepest  descent,  b)  Experimental  time  vs.  surface  tension  data  for  8  mMolar 
DLPC  with  exponential  fit  as  predicted  by  Equation  13b.  The  exponential 
behavior  of  t(g)  through  t  =  120  minutes  (dashed  line)  exhibits  remarkably 
good  agreement  with  experimental  data.  After  2  hrs  the  rate  of  surface  tension 
decay  slows  dramatically  and  approaches  the  asymptotic  limit  characteristic  of 
tightly  packed  monolayers.  See  text  for  details. 

Figure  8.  Dynamic  surface  tension  data  for  DLPC,  DMPC,  DPPC,  and  DSPC  under 
ambient  conditions  (22°  C)  and  at  a  bulk  concentration  of  1  (iMolar.  Also 
appearing  are  the  assymptotic  equilibrium  surface  tensions  of  these  systems. 
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Figure  9.  Dynamic  surface  tension  data  of  DLPC  in  the  low  concentration  limit  at  several 
different  temperatures.  Also  plotted  are  the  exponential  decays  predicted  by 
Equation  12  Discontinuities  in  the  38°  C  data  result  from  the  gradual  growth  and 
sudden  release  of  carbon  tetrachloride  vapor  bubbles.  From  the  exponential  fits 
come  rate  constants  for  monolayer  formation  according  to  Eq.  14. 

Figure  10.  Arrhenius  plots  of  the  monolayer  formation  rate  constants  for  the  four 

phosphatidylcholines  examined  in  this  work.  Activation  energies  are  calculated 
from  straight  line  fits  to  the  data.  Note  the  discontinuity  in  the  DMPC  data  ^ 
where  the  slope  changes  at  Tc  for  the  DMPC  gel-liquid  crystalline  bilayer  melting 
transition.  See  text  for  details. 

Figure  11.  Proposed  energetics  of  the  monolayer  formation  process  at  an  aqueous:CCl4 

interface  from  a  solution  of  PC  vesicles.  The  ordinate  plots  the  internal  energies 
of  gel  and  liquid  crystalline  phase  vesicles  as  well  as  the  internal  energies  of  the 
different  stages  of  monolayer  formation.  Internal  energies  of  the  vesicles  and 
monolayers  are  not  quantitatively  related  on  this  scale.  Along  the  abscissa  lies 
the  vesicle  rupture/monolayer  formation  coordinate. 
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